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Abstract. We present new global fits of the constrained Minimal Supersymmetric Standard 
CN | Model (cMSSM), including LHC 1/fb integrated luminosity SUSY exclusion limits, recent 

£^ ■ LHC 5/fb constraints on the mass of the Higgs boson and XENON100 direct detection data. 

Q\ . Our analysis fully takes into account astrophysical and hadronic uncertainties that enter the 

analysis when translating direct detection limits into constraints on the cMSSM parameter 
space. We provide results for both a Bayesian and a Frequentist statistical analysis. We find 
that LHC 2011 constraints in combination with XENON 100 data can rule out a significant 
y—{ \ portion of the cMSSM parameter space. Our results further emphasise the complementarity 

of collider experiments and direct detection searches in constraining extensions of Standard 

■ Model physics. The LHC 2011 exclusion limit strongly impacts on low-mass regions of 
cMSSM parameter space, such as the stau co-annihilation region, while direct detection 

5^ " data can rule out regions of high SUSY masses, such as the Focus-Point region, which is 

unreachable for the LHC in the near future. We show that, in addition to XENON100 
data, the experimental constraint on the anomalous magnetic moment of the muon plays 
a dominant role in disfavouring large scalar and gaugino masses. We find that, should the 
LHC 2011 excess hinting towards a Higgs boson at 126 GeV be confirmed, currently favoured 
regions of the cMSSM parameter space will be robustly ruled out from both a Bayesian and 
a profile likelihood statistical perspective. 
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1 Introduction 

The CMS and ATLAS collaborations have recently published new results for the search for 
Supersymmetry (SUSY) [1-3], based on proton-proton collisions with a center-of-mass energy 
t/s = 7 TeV, and an integrated luminosity of 1 fb^ 1 , recorded during 2011. No significant 
excess above the Standard Model (SM) predictions was measured, so that new constraints 
on SUSY could be derived, in particular for what concerns the cMSSM (constrained Minimal 
Supersymmetric Standard Model, see Refs. [4-8] and references therein), a simple construc- 
tion of minimal SUSY widely studied in the literature, which has become the de facto default 
model to evaluate SUSY searches and, often, prospects for dark matter searches [9-12]. Addi- 
tionally, an update on the search for the SM Higgs boson based on 5 fb" 1 of data has recently 
been announced by the ATLAS and CMS experiments [13]. In particular, the ATLAS ex- 
periment has reported a new 95% exclusion region constraining the allowed mass range for 
the SM Higgs boson to 115.5 — 131.0 GeV, as well as a 2.3<r excess at a mass of 126 GeV. 

In Ref. [14], we presented global fits of the cMSSM including data from precision tests of 
the Standard Model (SM), constraints from cosmology on the dark matter relic abundance, 
LEP bounds, and constraints on SUSY from LHC 2010 data. Additionally we studied the 
impact of recent exclusion limit on the WIMP-proton spin-independent scattering cross- 
section < J ^_ p from the XENON100 experiment [15], a two-phase time projection chamber 
with a liquid xenon target, on the cMSSM parameter space. Similar studies of the impact of 
LHC and XENON100 data on global fits of the cMSSM and other SUSY models are found in 
Refs. [16-20]. Our approach differs from other studies, because we provide results for both 
a Bayesian and frequentist statistical analysis, and included for the first time the relevant 
uncertainties (astrophysical and related to hadronic matrix elements) for the analysis of direct 
detection data. 

In this work we update the global fits of the cMSSM presented in Ref. [14] to include 
LHC 1 fb" 1 exclusion limits and LHC 5 fb _1 constraints on the Higgs boson (See also 
Refs. [17, 18, 20] for a discussion of the impact of LHC 2011 SUSY searches on the cMSSM, 
as well as on more general SUSY models. See Ref. [21] for a study of the impact of a 
hypothetical Higgs discovery at 125 GeV on the cMSSM parameter space). In our analysis, 
we apply all the experimental bounds used in Ref. [14], most notably constraints on cr?o_ 
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derived from the XENON100 experiment. We use nested sampling as a scanning algorithm 
[22] and study both the resulting Bayesian posterior distributions and the frequentist best fit 
regions. The methodology used has been described in detail in previous papers [14, 23-25], in 
particular we closely follow the procedure presented in [14]. We do not repeat the description 
of the methodology here and instead focus on the discussion of new physical results from the 
tighter LHC constraints. For further details on the statistical and theoretical background of 
our approach we refer the reader to Ref. [14] and references therein. 

Constraints on the WIMP mass and spin-independent cross-section derived from di- 
rect detection data strongly depend on the parameters describing the dark matter halo 
model, namely, the WIMP velocity distribution and the local density [26]. Additionally, 
constraints derived on the cMSSM parameters from the XENON100 exclusion limit depend 
on the hadronic matrix elements, which parameterise the contributions of light quarks to the 
proton composition. None of these quantities is precisely known, and neglecting the uncer- 
tainty on these quantities in the modelling can lead to incorrect inferences for the cMSSM 
parameters. Our inference procedure consistently takes into account our (lack of) knowledge 
of these parameters. We account for both astrophysical and nuclear physics uncertainties en- 
tering in the calculation of the event rate in direct detection experiments, by treating them 
as nuisance parameters in the analysis, which in the end are marginalised/maximised over 
to correctly infer the cMSSM parameters. 

This paper is organised as follows. In Section 2 we present the theoretical and statistical 
framework for the analysis. In Section 3 we present the results, Section 4 is devoted to 
the discussion of the prospects for discovering SUSY in the cMSSM in different detection 
channels. We present our conclusions in Section 5. 

2 Theoretical and statistical framework 

The cMSSM is described by four continuous parameters: the universal mass terms for scalars 
mo and gauginos mj/2, the universal scalar trilinear coupling Aq and the ratio of the Higgs 
vacuum expectation values tan/3. There is an additional discrete parameter, the sign of the 
Higgsino mass parameter, which is fixed to sgn(/i) = +1 in the following, since this is the 
value favoured by current measurements of the anomalous magnetic moment of the muon [27]. 
Despite its small number of parameters the cMSSM contains several of the main phenomeno- 
logical features of SUSY. Additionally the Lightest Supersymmetric Particle (LSP) is usually 
the lightest neutralino x?, a superposition of the neutral gauginos and higgsinos, which is an 
excellent cold dark matter candidate. This is a consequence of the structure of the Renor- 
malization Group Equations, which tend to make the Bino mass parameter small at the 
electroweak scale. Similarly, the stau is the lighter of all the sleptons, and in fact, it can 
become the LSP if the universal scalar mass is small with respect to the gaugino mass. In 
the regions where the stau is only slightly heavier than the neutralino a coannihilation effect 
in the early Universe helps reducing the neutralino relic density, making it easier to reproduce 
the observable dark matter abundance. Finally, for large values of the trilinear parameter 
the lighter stop can also become the LSP or provide a similar coannihilation effect to the 
relic density of neutralinos [28]. 

Residual uncertainties on the values of some SM parameters can strongly influence 
the constraints derived on the model parameters and observables [29]. Therefore, four SM 
nuisance parameters (namely, the top and bottom masses and the electroweak and strong 
coupling constants) are included in the scan. Additionally, when including constraints from 
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direct detection searches in the analysis, large uncertainties enter via astrophysical and nu- 
clear physics quantities. In Ref. [14] we provided a detailed discussion on the origin of these 
uncertainties and their impact on the calculation of the direct detection recoil rate. In the 
following we simply recall the resulting parameterisation of the uncertainties and refer the 
reader to the above paper for further details. 

The main astrophysical uncertainties are quantified using four nuisance parameters, 
describing both the abundance of dark matter and the WIMP velocity distribution. We use 
a simple Ansatz to describe the velocity distribution function (see section 3.2 in Ref. [14]). 
This parameterisation can correctly reproduce the different WIMP halo profiles that result 
from N-body simulations [30]. However, it is important to keep in mind that the actual 
local dark matter velocity distribution function remains unknown. N-body simulations have 
shown that the velocity distribution strongly depends on the detailed assembly history of 
our galaxy [31]. Additionally, for simplicity we neglect the presence of a dark disk [32-34] in 
our modelling of uncertainties. We do not include microhalos and localised velocity streams 
in our analysis, since the impact of these on direct detection signals has been estimated to 
be small (see e.g. Ref. [31, 35]). For a discussion of the different sources of astrophysical 
uncertainties affecting direct detection experiments see e.g. Ref. [36]. 

The hadronic uncertainties arise when translating constraints on the cMSSM parame- 
ters into constraints on the WIMP-proton spin-independent cross-section °"^o_ p - This cross- 
section depends on the hadronic matrix elements fx q = rn p < p\m q qq\p >, which represent 
the contributions of the light quarks q = {u, d, s} to the proton composition. These quantities 
are not precisely known and are treated as nuisance parameters in the analysis. Therefore, 
in addition to the four cMSSM model parameters, we also scan over up to eleven nuisance 
parameters, thereby fully accounting for those astrophysical and hadronic uncertainties. The 
specific ranges used for all the above nuisance parameters can be found in Table II of Ref. [14] . 

We investigate the constraints that can be placed on the cMSSM parameters using 
Bayesian methods to explore both the posterior probability density function (pdf) and the 
profile likelihood [25]. Bayesian inference is based on Bayes' theorem [37] 

which states that, given a likelihood function p(D|0) for a set of experimental data D 
and parameters of interest 0, and some prior pdf p(&) representing any information about 
the parameters available before taking into account the data, one can find the posterior 
probability p(0|D) for the parameters. For the purpose of parameter inference the Bayesian 
evidence p(D) acts as a normalisation constant, and is not important in the following analysis. 
As can be seen from Eq. (2.1) the posterior distribution depends on both the likelihood 
function and the prior. In the ideal case the experimental data are constraining enough to 
overcome the effect of the prior, such that the posterior distribution will be dominated by 
the likelihood function. When applying Bayes' theorem to complicated multi-dimensional 
models the data is often not sufficiently constraining to completely overcome the priors. By 
repeating the analysis for different choices of prior distributions one can assess to which 
extent the posterior pdf is dominated by the experimental constraints. Should the posterior 
display some residual prior dependence one has to be careful with the interpretation of the 
resulting constraints on the observables [24, 38]. Such a prior dependence has for example 
been observed in [14], where two different sets of priors were applied. We thus repeat each of 
our scans for two different sets of priors: "flat" priors (uniform in the cMSSM masses), and 



- 3 - 



"log" priors (uniform in the log of the cMSSM masses). Both sets of priors are uniform on Aq 
and tan (5. For the nuisance parameters we choose informative Gaussian priors, with the mean 
and standard deviation chosen appropriately to reflect current experimental constraints. The 
range covered by the scan for the cMSSM parameters, and the mean and standard deviation 
for the Gaussian priors on the nuisance parameters are given in Table I and Table II of 
Ref. [14], respectively. For scans that do not include astrophysical and hadronic uncertainties 
the nuisance parameters were fixed to their mean value. 

In order to identify constraints derived on a subset of one or two parameters one can 
consider two different quantities. The Bayesian marginalised one-dimensional posterior pdf 
for some parameter ©j can be found from the full posterior pdf by integrating out the 
remaining parameters 

p(0i|D) = J p(e|D)dei...de i _ide i+1 e n , (2.2) 

and equivalently for the 2-dimensional marginalised posterior pdf. Alternatively one can 
consider the frequentist profile likelihood function. The one-dimensional profile likelihood 
function for some parameter 0j can be found by maximising over all other parameters 

C{Qi) = max £(©). (2.3) 

fe>l,...,fe>;_l,te>;+i,...,fe) n 

These two statistical quantities have a different meaning and may lead to different conclusions. 
The profile likelihood function is ideal to discover small regions in parameter space that 
correspond to a large likelihood value. In contrast, the marginalised posterior pdf accounts 
for volume effects by integrating over hidden dimensions, and thus peaks around the region 
of highest posterior mass. Both of these quantities provide valuable information about the 
parameter space of interest, therefore we present results for both the marginalised posterior 
pdf and the profile likelihood function. To obtain the posterior distribution we modified 
the publicly available SuperBayeS vl.5.1 package [40] . As a scanning algorithm we use 
MultiNest v2.8 [22, 51], the running parameters are the same as in Ref. [14]. Each scan 
is repeated for both log and flat priors. The profile likelihood, which is in principle prior 
independent, is derived from combined chains of these two scans, as advocated in Ref. [25]. The 
chains used to derive the constraints on the cMSSM parameters and the observables are 
generated from approximately 7 (43) million likelihood evaluations, with an overall efficiency 
of 0.1 (0.02) for fixed (varying) astrophysical and nuclear physics nuisance parameters. 

The experimental data included in the likelihood function are given in Table III of Ref. 
[14]. The full likelihood function includes precision tests of the SM, constraints on the dark 
matter relic density from 7-year WMAP data and constraints from accelerator searches. We 
also include recent data from the XENON100 direct detection experiment, obtained with an 
effective volume of 48 kg and an exposure of 100.9 days between January and June 2010 [15]. 
Three candidate scattering events were detected within the signal region, which is compatible 
with the expected number of background events b = 1.8 ± 0.6, resulting in tight exclusion 
limits in the (rn^o.azk ) plane. Our full likelihood function is described in detail in Ref. 
[14], which we refer the reader to for further details. 



For this paper, the public SuperBayeS vl.5.1 code has been modified to interface with with DarkSUSY 
5.0.5 [41, 42], SoftSUSY 2.0.18 [43, 44], MicrOMEGAs 2.0 [45, 46], Superlso 2.4 [47, 48] and SusyBSG 1.4 
[49, 50]. 
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In this paper we upgrade the constraints from SUSY searches at the LHC to include 
recent exclusion limits from 1 fb _1 integrated luminosity data, as well as the latest limits 
on the SM Higgs mass obtained with 5 fb _1 integrated luminosity worth of data. We focus 
on constraints presented in Ref. [1] by the CMS collaboration, which currently provides the 
most stringent constraints on SUSY. The CMS data set was collected for a centre-of-mass 
energy of y/s = 7 GeV and an integrated luminosity of 1 fb _1 . A SUSY signal is searched for 
in hadronic events with two or more jets and missing transverse energy using the kinematic 
variable olt- No significant signal beyond the SM predictions was observed, and so lower 
limits on the cMSSM masses mo, mi/2 were derived, as shown in Fig. 5 of Ref. [1]. While 
this exclusion curve was obtained for fixed values tan /3 = 10, Aq = 0, the decay channels used 
to obtain this limit are relatively insensitive to the values of these two parameters. Thus, 
we can treat the LHC constraint as approximately independent of tan f3 and Aq [52] . The 
LHC 2011 exclusion limit is included in the likelihood function by defining the likelihood of 
samples corresponding to masses below the limit to be zero. 

The ATLAS and CMS collaborations have recently announced new stringent constraints 
on the Standard Model Higgs boson mass, 115.5 < m/j < 131 GeV (at 95%, from a combined 
11 channels analysis of ATLAS data), derived from searches with 5 fb _1 integrated luminosity. 
Furthermore, a ~ 3.6a local excess has been reported for ~ 126 GeV; the significance 
is reduced to ~ 2.3cr when the "look-elsewhere" effect is considered. Although it is too 
early to interpret such an excess as an indication of the Higgs, we include the updated 
lower limit in our analysis, by applying the 95% lower limit of m/j > 115.5, with an added 
theoretical uncertainty of 3 GeV, according to the procedure described in Ref. [23]. Although 
the ATLAS limit applies to the Standard Model Higgs, the lightest Higgs in the cMSSM is 
almost invariably SM-like, hence we can adopt the above limit in this context. 

3 Results 

The results of our numerical analysis will be presented in three steps. First we show the im- 
pact on the cMSSM parameter space resulting from all the constraints described in the previ- 
ous section (except those from direct detection experiments), including the latest LHC data. 
Afterwards we combine these results with XENON100 data, fully marginalising/maximising 
over the astrophysical and hadronic nuisance parameters. Finally, we discuss in some detail 
the impact of the experimental constraint on the anomalous magnetic moment of the muon 
on the cMSSM parameter space. 

3.1 Impact of LHC 2011 data 

In Fig. 1 we show the impact of LHC 2011 data on the cMSSM in the {jni/2i m o) plane 
(left), the (tan/3,Ao) plane (centre) and the (m^o,a%_ ) plane (right). The two-dimensional 
marginalised posterior pdfs are given for both flat (top panel) and log (central panel) priors. 
The profile likelihood is shown in the bottom panel. Blue contours correspond to constraints 
derived without including LHC data, black contours show constraints on the parameters after 
the latest LHC data are included. 

For the pre-LHC contours we observe several regions in the (mo, ^1/2) plane that are of 
particular interest. The first is the so-called focus point (FP) region [54, 72], or hyperbolic 
branch [55], which corresponds to a region in parameter space where the Higgsino mass 
parameter fi is very small, while scalar masses can be very large. Due to the smallness of \i 
the neutralino has a large higgsino component, which leads to its relic density to decrease 
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Figure 1. Black contours: posterior pdf (upper panels, for flat and log priors) and profile likelihood (lower 
panels) for the cMSSM parameters, including all present-day constraints (WMAP 7-years, LHC 1 fb _1 SUSY 
searches and 5 fb -1 Higgs limits included), except XENON100. From the inside out, contours enclose 68%, 
95% and 99% of marginal posterior probability (top two rows) and the corresponding profiled confidence 
intervals (bottom panels). The black cross represents the best fit point, the black dot the posterior mean 
(for the pdf plots). Parameters describing astrophysical and hadronic uncertainties have been fixed to their 
fiducial values. Blue contours represent the constraints obtained without the inclusion of LHC data. In the 
plots on the left, the dashed/green line represents the current LHC exclusion limit, while in the right-most 
plots the red/dashed line is the 90% exclusion limit from XENON100, from Ref. [15], rescaled to our fiducial 
local DM density of p loc = 0.4 GeV/cm 3 . 

[56], such that current cosmological constraints on the dark matter relic abundance can be 
satisfied. In the (rrii/2, n^o) plane it corresponds to a large area at sizable rriQ > 1 TeV and 
relatively small nii/2- 

A second region of interest is the stau co-annihilation (SC) region. In this region the 
lightest stau is slightly heavier than the neutralino LSP. As a result, the neutralino relic 
density is reduced by neutralino-stau co- annihilations, so that the WMAP constraint can be 
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LHC 2011 


LHC 2011 + XENON100 


LHC 2011 w/o Sa% uay 


Input parameters 


m [GeV] 


282.19 


267.54 


188.80 


%/2 [GeV] 


691.76 


635.56 


908.06 


[GeV] 


685.35 


935.04 


-630.62 


tan j3 


33.74 


29.75 


8.54 


M t [GeV] 


174.108 


173.397 


173.338 


m b {m b ) Kls [GeV] 


4.214 


4.221 


4.217 




127.952 


127.952 


127.955 


a s (M z ) MS 


0.116 


0.119 


0.117 


Observables 


m h [GeV] 


117.1 


115.6 


119.4 


m^o [GeV] 
6a^ SY x 10 10 


287.2 
22.82 


262.2 
32.57 


382.5 
2.32 


5/2(5 -> V s7 ) x 10 4 


3.04 


3.42 


3.36 


A _ x 10 2 


8.42 


8.36 


8.04 


Bi2(D fl -> ri/) x 10 2 


4.81 


4.81 


4.82 


Bi2(Z? fl -»• /ii/) x 10 3 


4.97 


4.97 


4.98 


BR(D -> /xi/) x 10 4 


3.86 


3.86 


3.86 


°§-p [pb] 

^ s !? [ P b] 


2.6 x 10~ 9 


3.8 x 10" 9 


3.0 x 10~ 10 


3.8 x 10- 7 


5.8 x 10~ 7 


3.7 x 10~ 8 


^o/l 2 


0.1174 


0.1121 


0.1118 



Table 1. Input cMSSM parameters (and some other interesting quantities) for the best fit points. 
The column "LHC 2011" is for the case where LHC 2011 data, but no direct detection data is applied, 
while "LHC + XENON100" denotes the case where XENON100 data is added to the analysis and 
astrophysical and hadronic nuisance parameters have been included in the scan and profiled over. 
The column "LHC 2011 w/o 5a^ USYv is for the case where all constraints, including LHC 2011 data, 
except for direct detection data and the Sa^ u constraint, are applied. 

satisfied. For fixed Aq and tan/3 the SC region appears as a narrow band on the (m^2i"io) 
plane spanning a range of gaugino masses m^ 2 ~ a few hundred GeV at small scalar masses. 
For varying Aq and tan f3 this region appears slightly smeared. 

A third possibility for neutralinos to reproduce the correct relic abundance is for them 
to undergo a resonant annihilation mediated by a relatively light pseudoscalar Higgs. This is 
possible when m^o ~ Im^o and typically happens for large values of tan f3 and relatively large 
values of both the gaugino and scalar mass parameters. The position of this so-called Higgs 
funnel region in the (min, tuq) plane strongly depends on other parameters (in particular on 
tan j3 and the trilinear term) . 

Finally, there is a narrow viable region in which the lightest neutralino mass is approx- 
imately half of the lightest Higgs mass, m^o ~ 60 GeV. In this case, the s-channel Higgs 
exchange in the neutralino annihilation cross section in the early Universe is resonant and 
the resulting relic abundance is significantly reduced. This appears in Fig. 1 as a narrow 
area with ~ 100 GeV and spanning several orders of magnitide in tjiq. For a further 
discussion of the features of the cMSSM see e.g. Ref. [24]. 
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LHC 2011 


LHC 2011 + XENON100 


LHC 2011 w/o 6a?. ulSY 


Observable 








SM nuisance 


1.059 


0.528 


0.100 


Astro nuisance 


N/A 


1.052 


N/A 


Hadronic nuisance 


N/A 


0.172 


N/A 


M w 


0.751 


1.074 


1.228 


sin 2 e ff 


0.158 


0.084 


0.044 




O.DbO 


0.12 1 


0.000 


BR(B -> X sl ) 


1.665 


0.113 


0.218 


AM Bs 


0.331 


0.243 


0.306 


J3 iX ^ iJ u Tils J 

BR(B u ^rv) SM 


0.918 


0.867 


0.551 








A - 


3.305 


3.232 


2.812 


BR(B~ j >Dtu) 
BR(B^-Deu) 


0.868 


0.859 


0.799 




0.390 


0.382 


0.328 


BR(D S -> ri/) 


2.252 


2.248 


2.214 




3.136 


3.131 


3.098 


£i?(D -»• pi/) 


0.008 


0.008 


0.008 


ft^o/l 2 


0.171 


0.0002 


0.002 




1.07 


2.157 


0.287 


XENON100 


N/A 


0.019 


N/A 


LHC 


0.0 


0.0 


0.0 


Sparticles (LEP) 


0.0 


0.0 


0.0 


BR(B S -»• 


0.0 


0.0 


0.0 


Total 


16.74 


16.29 


12.00 


X 2 /dof 


1.86 


1.81 


1.50 



Table 2. Breakdown of the total x 2 by observable for the best fit points, for both the case where 
all data, including LHC 2011 constraints, except for XENON100 data were applied (left column), 
when all data including XENON100 results were applied and nuisance parameters were profiled over 
(central column), and when all data including LHC 2011 constraints except for XENON100 data and 
the Sa^ u constraint were applied (right column). 

The LHC 2011 data has a strong impact on the posterior distribution in the (miw, mo) 
plane, as shown on the left-hand side of the upper and central panels in Fig. 1. The LHC 
2011 exclusion limit cuts deep into the stau co-annihilation region. A large region that 
was previosuly favoured at the 68% and 95% level (see Ref. [14]) is disfavoured with high 
confidence and the posterior contours are pushed towards higher values of tuq and 
Notice that the narrow band with resonant Higgs annihilation remains viable for neutralinos 
with a mass around 60 GeV. 

As can be seen on the right-hand side of Fig. 1 the LHC exclusion limit also has impor- 
tant consequences for direct detection of the cMSSM: The exclusion of small gaugino masses 
leads to small m^o being strongly disfavoured, such that a large part of previously favoured 
cMSSM parameter space at small and intermediate WIMP masses is now excluded at the 
99% level. In the (mjo,a?o_ ) plane the FP region appears as a large island corresponding 

to relatively large a%_ , lying above the projected XENON100 90% exclusion limit shown 
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in red (not applied here). One can already see from this limit that adding direct detection 
data has the capability of ruling out significant portions of the FP region, which appears at 
high scalar masses and is therefore unaffected by current LHC constraints. 

In the central plot one can see that LHC 2011 constraints have a minimal impact on 
the credible regions for tan/3 and Aq. This is to be expected because, as stated above, the 
LHC exclusion limit is fairly independent of the precise values of these two parameters. 

These results can be compared to the impact of earlier LHC 35 pb _1 limits on global fits 
of the cMSSM, as presented in Ref. [14]. Qualitatively the resulting posterior distributions are 
quite similar. The main difference is that while the 35 pb" 1 exclusion limit had a fairly modest 
impact on the cMSSM parameter space, the LHC 2011 data leads to exclusion of higher values 
of the scalar and gaugino masses, so that a large portion of the stau-coannihilation region 
is ruled out. LHC 2011 data also impacts on the posterior distribution in the (m^o,cr|o_ p ) 

plane, which was more or less unaffected by the 35 pb^ 1 data. Now a significant range of 
small WIMP masses is excluded at low and intermediate a. Q 



The above observations are true for both the log and the flat prior case. Due to "volume 
effects" the flat prior contours extend to larger scalar and gaugino masses, otherwise the 
posterior pdfs resulting from the two scans agree well in all three planes. This is encouraging, 
since it indicates that experimental data is becoming constraining enough to dominate the 
inference results, which are now only weakly dependent on the choice of prior, although some 
residual volume effects remain in the relative posterior weight of the FP region with respect 
to other parts of the parameter space. 

The 2D profile likelihood results 2 are shown in the bottom panels of Fig. 1. Similarly to 
the Bayesian analysis, the LHC 2011 exclusion limit cuts deep into the stau-coannihilation 
region. The favoured region of parameter space is still located at relatively small scalar and 
gaugino masses, not far above the new LHC exclusion limit. The profile likelihood results 
are qualitatively very similar to the Bayesian results, most importantly the FP region is 
still allowed at 99% confidence level. The 99% contours for the profile likelihood are very 
similar to the extent of the 99% Bayesian credible intervals shown in the upper panels. As 
we already pointed out in Ref. [14], the 99% confidence region is much larger than expected 
from the 68% and 95% contours, indicating that the tails of the profile likelihood function 
are highly non-Gaussian. This demonstrates that the use of a high resolution scan such as 
the one we employ in this paper is required in order to correctly map out the shape of the 
profile likelihood function and accurately infer the extent of the 99% confidence level. 

We now examine the best fit points in more detail. Details of the best fit points are 
given in Table 1, a breakdown of the total x 2 by observable is given in Table 2, where \ 2 can 
be found by summing over the likelihood Li for each observable 



For observables that are constrained by exclusion limits we choose In Ci = for points outside 
the allowed region. An exception is the XENON100 exclusion limit, as described in detail 

2 In order to construct the profile likelihood contours presented here, it was necessary to remove from our 
scans the highest likelihood point, which showed an anomalously low x 2 . This could be traced back to the 
fact that this point had a very small mass difference between the second-heaviest neutralino and the heaviest 
smuon, with a fractional difference of less than 10 -6 . This led to reproducing very closely the observed 
anomalous magnetic moment of the muon, but with a large amount of numerical fine tuning. We removed 
this finely-tuned point from the scan. 



X°-P 
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in Ref. [14]. For Gaussian-distributed observables we normalize the likelihood such that 
ln£j = when the theoretical and the experimental values match. In this analysis we use 
13 "active" Gaussian data points (see Table 3 in Ref. [14]) and scan over 4 free cMSSM 
parameters. Therefore, in total there are 9 degrees of freedom. 

The best fit point for this scan is found in the SC region, corresponding to small scalar 
and gaugino masses < 1 TeV that are still allowed by LHC 2011 constraints. The observables 
having the largest contribution to the total x 2 m Table 2 are the isospin asymmetry Ao- and 
the branching ratios of the BR{D S — > tu) and BR{D S — > nv) processes which are sensitive 
to new physics mainly through the mass of the charged Higgs bosons and to some extend 
to tan/3 [57], in particular at low mjj±. Our best fit points have relatively large mjj± and 
therefore the theoretical values for those processes are SM-like and have a discrepancy at 
the ~ 1.5cj level with the experimental values. For the isospin asymmetry Ao- new physics 
contributions become important mainly to low mi/2 values [58]. Hence the predictions for our 
best fit points are again SM-like, showing a discrepancy at the 2 a level with the experimental 
central values. 

3.2 Impact of XENON100 data, including hadronic and astrophysical uncer- 
tainties 

The effect of adding the XENON100 exclusion limit is shown in Fig. 2. Astrophysical and 
hadronic uncertainties are accounted for by marginalising/profiling over the corresponding 
nuisance parameters. We have checked that, in agreement with the findings of Ref. [14], the 
overall impact of including these uncertainties in the analysis leaves our results qualitatively 
unchanged, though quantitatively it affects the overall limits. 

As can be seen in the upper two panels of Fig. 2, the inclusion of XENON100 data 
has a significant impact on the two-dimensional posterior distributions for the cMSSM. The 
focus point region, which could previously not be excluded at the 99% level, is now strongly 
disfavoured; the parameter space is significantly reduced. Aside from residual volume effects 
the contours in the (m^.o,o"?o_ p ) plane for the two different choices of priors are in relatively 
good agreement with each other. Note that for the posterior with flat priors the SC region 
is strongly disfavoured. This is not a physical effect, as can be seen from the fact that 
the combined best fit point from the flat and log prior scans is found inside the SC region. 
Instead, this is a result of an underexploration of the cMSSM parameter space by the flat 
prior scan. Due to the high dimensionality of the parameter space small regions of high 
likelihood can easily be missed. The flat prior scan is particularly vulnerable to this, since it 
explores the low mass regions in much less detail than the log prior scan. 

By comparing the top and central right-hand panels in Fig. 2 we can see that the 
99% contours are very similar for the posterior pdfs for both choices of priors. Thus, the 
constraints on m^o and a ^_ p are becoming increasingly independent of the choice of priors. 

The PL results, shown in the bottom panels, qualitatively agree well with the Bayesian 
analysis. The parameter space shrinked significantly. The contours are tighter than the 
Bayesian credible regions, which is explained by the high likelihood value of the best fit 
point (see table 2). This illustrates that profile likelihood results have to be interpreted with 
care; the extent of the contours strongly depends on the likelihood value of a single point in 
parameter space found by the scan. Again the best fit point is found in the SC region. Note 
that this best fit point corresponds to a relatively low rrih, such that a sizeable contribution 
to the x 2 arises from the 5 fb" 1 95% Higgs exclusion limit. Similar to the results of the 
Bayesian analysis, the FP region is ruled out at 99% confidence level by the XENON100 
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Figure 2. As in Fig. 1, but now black contours include XENON100 data (considering hadronic and 
astrophysical uncertainties as nuisance parameters), while the blue empty contours show for comparison the 
case where no direct detection data are included (from the inside out: 68%, 95% and 99% regions). We observe 
a strong suppression of the viability of the FP region. Notice that the XENON 100 90% limit (red/dashed 
line) has been included only to guide the eye, as our implementation of the XENON100 data is slightly more 
conservative than the procedure adopted in Ref. [15]. 



data, as can be seen explicitly in the rightmost plot. This clearly illustrates the potential of 
direct detection experiments to constrain SUSY. 

3.3 Impact of the 5a^ u constraint 

The muon anomalous magnetic moment provides an interesting window to new Physics, 
since it is very accurately measured. A constraint on the supersymmetric contribution to 
this observable, 5a^ USY , can be extracted by comparing the experimental result [59], with 
the theoretical evaluations of the Standard Model contribution [60-62]. Although the latter 
have become increasingly precise in the last decade, they are still subject to theoretical uncer- 
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Figure 3. As in Fig. 1, but now black/filled contours result from scans that do not include the experimental 
constraint on the anomalous magnetic moment of the muon. Blue contours correspond to the results obtained 
when applying this constraint, and thus are identical to the black contours in Fig. 1. 

tainties, especially in the computation of the hadronic loop contributions. According to the 
most recent evaluations, when e + e~ data are used the experimental excess in = (g^ — 2) /2 
would constrain a possible supersymmetric contribution to be 5a^ USY = (29.6±8.1) x 10~ 10 , 
where theoretical and experimental errors have been combined in quadrature. However, when 
tau data are used a smaller discrepancy (2.4 a) with the experimental measurement is found 
[61]. 

In Ref. [24] it has been shown that the preference for small ttlq and m^ 2 in global 
fits of the cMSSM is strongly driven by the Sa^ USY constraint. In order to evaluate the 
dependence of our results on this constraint, we repeat the analysis presented in Section 3.1 
after dropping the experimental constraint on 5a^ USY . The results are shown in Fig. 3. For 
comparison, blue contours show the constraints derived on the cMSSM when including the 
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Figure 4. ID marginal pdf for flat priors (thin solid/blue), log priors (thick solid/red) and ID profile 
likelihood (dashed/black) for the gluino mass, the neutralino mass, the SI and the SD scattering cross section 
(from left to right). Top panels include all data except XENON100, bottom panels include XENON100 
data with hadronic and astrophysical uncertainties fully marginalised/maximised over. The best fit point is 
indicated by the encircled black cross. 

constraint in the analysis. 
As can be seen the posterior contours for both sets of priors are significantly expanded 
outwards towards larger values of tuq and m^j. The FP region is allowed at the 68% level. 
Clearly, the 5a^ USY constraint is the dominating factor for the exclusion of high scalar masses 
mo ~ a f* ew TeV. No other constraints strongly disfavour these regions. However, even when 
excluding the 8a^ USY constraint from the scan, the posterior probability for small scalar 
and gaugino masses remains high; these mass scales are favoured by a number of different 
constraints. 

These effects are even more pronounced for the profile likelihood analysis (bottom panels 
of Fig. 3). The extent of all three contours increases significantly when the 5a^ USY constraint 
is excluded from the scans. High scalar and gaugino masses are allowed at the 68% level. 
Without the constraint on the anomalous magnetic moment of the muon no conclusions can 
be drawn on the allowed ranges of tan/3 and Aq. The allowed region of parameter space in 
the (m^,o,(T?o_ p ) plane also increases significantly. The x 2 /dof for the best fit point, shown 

in Table 2, is significantly reduced, which emphasises that the 5a^ USY constraint is in strong 
conflict with several other constraints. However, the best fit point is still found at small ttiq 
and mi/2, such that the SC region remains the most favoured region even when dropping the 
constraint. 

Therefore, global fits of the cMSSM are still strongly driven by the Sa^ USY constraint, 
which significantly disfavours large scalar and gaugino masses and therefore pushes both the 
posterior and the profile likelihood contours to small values of uiq and m^. 



- 13 - 



4 Updated prospects for cMSSM discovery and Higgs excess 

We now present the implications of the above results for the prospects of detecting the 
cMSSM using the LHC, direct detection and indirect detection experiments. In Fig. 4 we 
show the ID marginal posterior distributions for both log (thick solid/red) and flat (thin 
solid/blue) priors and the ID profile likelihood functions (dashed/black) for the gluino mass 
nig , the mass of the lightest neutralino m^o , the spin-independent scattering cross section 
<7?o_ p and the spin-dependent neutralino-proton scattering cross section £7 |?_ p - 

The ID marginal posterior distributions for the log prior scan are quite similar to the 
ID profile likelihood results. The ID posterior pdf for the flat prior scan can differ from these 
distributions, especially for the gluino and the neutralino masses. This is due to the flat prior 
suffering from residual volume effects. Comparing these results with the corresponding plots 
for the LHC 2010 data presented in Ref. [14] it can be seen that now larger gluino masses 
rrig > 1500 GeV are favoured. The only exception to this is the very narrow area where 
60 GeV neutralinos have a resonant annihilation through s-channel diagrams mediated by 
the lightest Higgs. The gluino mass in those cases is around 450 GeV. This is apparent in 
the ID plots for the gluino and neutralino masses in Fig. 4. 

A sizable range of neutralino masses m x ~ 100 — 200 GeV is excluded when combining 
LHC 2011 and XENON100 constraints. This includes part of the mass region which was 
previously favoured by the LHC 2010 results. Additionally, cross-sections <r?o_ larger than 

Xi P 

10 -8 pb are strongly disfavoured. However, prospects for detection of the cMSSM by the 
next generation of direct detection experiments remain good, with our best fit point easily 
within reach of the future XENON1T experiment. This experiment is expected to probe the 
vast majority of the current 99% contours by 2015 and exclude cross-sections °"~o_p > 10~ 10 
pb for WIMP masses up to 300 GeV (see e.g. Ref. [63]). The expected 90% exclusion limit of 
the XENON1T experiment is indicated on the bottom right-hand panel in Fig. 2. The LHC 
2011 limit also tightens the constraints on the spin-dependent cross-section (T^$_ p , strongly 

disfavouring °jp_ p > 10~ 6 pb for both the posterior and the profile likelihood contours. This 
further reduces indirect detection prospects of the cMSSM and makes its detection by the 
IceCube neutrino observatory, which is able to probe <r??_ > 3 x 10~ 5 pb (see e.g. Ref. 

Xi P 

[64]), highly unlikely. 

The effect of LHC 2011 constraints on the mass of the lightest Higgs are presented 
in Fig. 5. From left to right we show constraints on rrih including LHC 2011 limits but 
no direct detection data, LHC 2011 constraints combined with XENON100 data and LHC 
2011 limits without XENON100 data and excluding the constraint on 5a'^ USY . Compared to 
LHC 2010 results, small Higgs masses are now strongly disfavoured, so that the probability 
distributions move to larger m^. This is a result of the new 95% lower limit > 115.5 
GeV. In the left and central plots the profile likelihood and the posterior pdf for log priors 
peak relatively closely to this new lower limit, while the flat prior scan favours slightly higher 
values of m^. When excluding the constraint on the magnetic moment of the muon from the 
analysis slightly higher masses > 120 GeV are favoured from both the Bayesian and the 
profile likelihood statistical perspective (see right-hand panel in Fig. 5). 

While we did not implement the excess at ~ 126 GeV reported by the ATLAS 
collaboration explicitly, it is still obvious from Fig. 5 that, should this excess be confirmed, 
the cMSSM will be strongly disfavoured. Both the frequentist profile likelihood and the 
Bayesian posterior probability approach zero at mh = 126 GeV, independent of the choice of 
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Figure 5. ID marginal pdf for flat priors (thin solid/blue), log priors (thick solid/red) and ID profile 
likelihood (dashed/black) for the lightest Higgs mass mu- The results come from the implementation of all 
experimental data, including LHC 2011 data, except for direct detection constraints (left), all data including 
XENON100 data with astrophysical and hadronic uncertainties fully marginalised/maximised over (centre) 
and all data except direct detection data and excluding the Sa^ USY constraint. The best fit point is indicated 
by the encircled black cross. 

prior. This conclusion remains true even if the Sa^ USY constraint is excluded from the scan 
(see right panel of Fig. 5). 

At this point we comment on how we expect the results of our analysis to change when 
using priors other than the non-informative fiat and log priors applied in this work. 

In the literature Bayesian studies of the cMSSM have been performed which attempt 
to incorporate the SUSY naturalness criterium. Namely, SUSY soft-masses should not be 
far from the experimental electroweak (EW) scale in order to avoid unnatural fine-tuning 
to obtain the correct size of the EW symmetry breaking. In some studies a penalisation of 
the fine-tuned regions has been implemented, e.g. by using a conveniently modified prior 
for the cMSSM parameters [66, 67]. On the other hand, in Ref. [68] it has been shown 
that the naturalness arguments arise from the Bayesian analysis itself, with no need of in- 
troducing "naturalness priors" . The key is when the experimental value of Mz is considered 
in the same way as other experimental information (usually Mz is fixed to its experimen- 
tal value and the Higgsino mass parameter \x is predicted from the EW symmetry breaking 
conditions). Marginalising over [i results in a factor in the Bayesian posterior, where 
c M = |cHnM|/<91n/i| is the conventional Bar bier i-Giudice measure of the degree of fine- 
tuning [70, 71] (for details on this derivation see Ref. [68, 69]). This precisely agrees with 
the "naturalness prior" which is introduced by hand in Ref. [66]. Thus, the presence of 
this fine-tuning parameter in the denominator penalises the regions of parameter space cor- 
responding to large fine-tuning. As a result the only region with large soft-masses that is 
not disfavoured is the FP region, in which naturalness is preserved [72]. Indeed, this region 
contains a large portion of the Bayesian posterior probability in the presence of the DM 
relic abundance constraint, especially when the constraint on the anomalous magnetic mo- 
ment of the muon is excluded from the analysis [69]. As was shown above, the addition of 
XENON100 data strongly disfavours the FP region, therefore one would expect the bulk of 
the posterior probability to fall within the low and intermediate soft-masses region, leading 
to similar conclusions as the ones resulting from our log prior scan. 
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5 Conclusion 

In this work we have presented new global fits of the cMSSM, including exclusion limits 
from LHC 2011 SUSY and Higgs searches and XENON100 direct detection data. In addition 
to uncertainties on Standard Model quantities, our analysis takes into account both astro- 
physical and hadronic uncertainties that enter the analysis when translating direct detection 
limits into constraints on the cMSSM parameter space. XENON100 data in combination 
with LHC 2011 constraints can rule out a significant portion of the cMSSM parameter space. 
In addition to the FP region being ruled out at the 99% level by XENON100 data, we have 
shown that the size of the SC region is significantly reduced by LHC 2011 constraints on mo 
and mi/2- LHC 2011 data strongly impacts on the cMSSM parameter space, pushing allowed 
particle masses to higher values and reducing the favoured values of both spin-dependent and 
spin-independent neutralino-proton cross-sections. 

This has important implications for direct detection of the cMSSM, with LHC con- 
straints cutting into the low m^o region of cMSSM parameter space and, in combination 

with direct detection constraints, pushing the contours to lower a??, . Our results further 

' r to x{-p 

emphasise the complementarity of collider experiments and direct detection searches. While 
LHC data are placing tight constraints on large portions of cMSSM parameter space, direct 
detection data is crucial to rule out regions of high SUSY masses, such as the FP region, 
which is unreachable for the LHC in the near future (see also the discussion in Ref. [65]). 

We further investigated the importance of the experimental constraint on the anomalous 
magnetic moment of the muon. We found that this constraint is the dominating reason for 
large scalar and gaugino masses being disfavoured in our analysis. All other constraints (aside 
from XENON100) leave this range of masses fairly unconstrained. 

Finally, we investigated constraints on the lightest Higgs mass resulting from LHC 2011 
and XENON100 data, is pushed towards larger values due to the new LHC 2011 95% 
confidence level. We have qualitatively shown that, should the excess around ~ 126 
GeV be confirmed by future data sets, the residual cMSSM parameter space will be strongly 
disfavoured. 

While this work was under completion, similar analyses of the constraints on the cMSSM 
parameter space resulting from LHC 2011 SUSY searches and XENON100 data were per- 
formed in Ref. [18, 20]. Our conclusions are qualitatively similar to Ref. [20], while in Ref. [18] 
it was claimed that the FP region can not be excluded at 95% level by XENON100 data and 
a negative outlook on direct detection prospects of the cMSSM was given. This is a result 
of a very conservative estimation of the XENON100 exclusion limit in Ref. [18]. Our anal- 
ysis differs from Ref. [20], where astrophysical uncertainties are not included in the analysis 
and results are presented for the profile likelihood only. In contrast, we provide a detailed 
comparison between Bayesian and profile likelihood results. Ref. [18] provides a similar com- 
parison, however, the scanning parameters chosen are not suitable for an accurate mapping of 
the profile likelihood function (according to the prescriptions given in Ref. [25]), so that the 
extent of the resulting confidence levels may vary strongly from scan to scan. Additionally, 
in Ref. [18] the impact of astrophysical and hadronic uncertainties is taken into account by 
smearing out the XENON100 exclusion limit, instead of marginalising/maximising over the 
corresponding nuisance parameters. In contrast to both Ref. [18] and Ref. [20] we include 
recent 5 fb" 1 LHC constraints on the Higgs mass, so that we provide a more sophisticated, 
up-to-date statistical analysis of the cMSSM parameter space. 
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